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Neutron scattering study of Kondo lattice antiferromagnet YbNiSia 
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The Kondo lattice antiferromagnet YbNiSi3 was investigated by neutron scattering. The 
magnetic structure of YbNiSi3 was determined by neutron diffraction on a single-crystalline 
sample. Inelastic scattering experiments were also performed on a pulverized sample to study 
the crystalline electric field (CEF) excitations. Two broad CEF excitations were observed, 
from which the CEF parameters were determined. The temperature dependence of the mag- 
netic susceptibility x an d the magnetic specific heat C m ag were calculated using the deter- 
mined CEF model, and compared with previous results. 
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1. Introduction 

Recently Yb-based compounds have been investigated intensively owing to their variety 
of unusual physical properties resulting from a competition between the Kondo effect and 
long-range inter-atomic exchange interactions, together with crystalline electric field (CEF) 
effects. The Kondo effect leads to a Fermi-liquid state at low temperatures made of strongly 
correlated heavy quasiparticles with enhanced effective masses, called heavy fermions. On the 
other hand, long-range inter-atomic exchange interactions tend to yield magnetically ordered 
states. In Yb-based heavy fermion systems, CEF effects split the Yb ion multiplet and play an 
important role in the physical properties. 1 To understand the behaviors of Yb and 4/ electron 
systems, more systematic and in-depth studies are imperative. 

We have studied YbNiSi3 using neutron scattering techniques. YbNiSi3 is an orthorhom- 
bic Kondo lattice heavy-fermion compound, which shows an antiferromagnetic state below 
Tj\j = 5.1 K. 2 This transition temperature Tn = 5.1 K is especially high among Yb 
compounds, indicating a strong magnetic exchange interaction. YbNiSi3 is thus considered 
a model case of almost +3 valence Yb-based heavy-fermion antiferromagnetic compounds, 
which can help a full recognition of a competing relation between the magnetic interaction 
and the Kondo effect. 
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The effects of CEF excitation on the magnetic susceptibility \ was pointed by Avila et al? 
The magnetic moments were considered to align along the easy magnetic axis b from mag- 
netic measurements, and the effective magnetic moment was estimated to be smaller than 
2 m the ordered state. 2 To advance the understanding of YbNiSi3, we performed neu- 
tron scattering experiments in this present work. Elastic-scattering measurements using a 
single crystalline sample were performed to resolve the exact magnetic structure. In addition, 
inelastic-scattering measurements using a pulverized sample were carried out to observe the 
CEF excitation scheme. The CEF model was then applied to analyze and explain some macro- 
scopic physical properties which have already been reported. 2 ' 3 We also discuss the ground 
state of Yb ion with the CEF model analysis, which the multiplet of the ground state has 
been a subject of a considerable discussion. 2 ' 3 

2. Experiment 

2.1 Magnetic Structure 

Neutron diffraction measurements were performed on the triple-axis spectrometer PONTA 
installed by ISSP in the JRR-3 reactor of the Japan Atomic Energy Agency ( J AEA). The (002) 
reflection of pyrolytic graphite (PG) was used to monochromate and analyze the neutron 
beam, together with a sapphire and two PG filters to eliminate higher order contaminations. 
The neutron beam energy was fixed at 14.7 meV. The horizontal collimation of 40' - 40' - 
40' - 80' was employed. The single-crystalline sample, sealed in an Al can with He gas, was 
mounted in an ILL-type orange cryostat. The sample size was about 3 x 3.5 x 0.5 mm 3 . 

An uj — 29 scan was performed after moving to the peak profile center-of-mass of an u 
scan on each reciprocal lattice point studied, below and above the Neel temperature. The 
scan profiles were fitted with a gaussian function to estimate the integrated intensities. The 
integrated intensities were calculated by multiplying the full-width at half-maximum (FWHM) 
by the peak height of each reflection. The diffraction experiments were first performed on the 
scattering plane (0 k I) and then changed to the (h k 0) plane. 

2.2 Crystalline Electric Field 

Powder neutron inelastic-scattering measurements were also performed on PONTA oper- 
ated in the neutron energy-loss configuration. The powder sample was obtained by pulverizing 
about 1.1 grams of single crystals. The final neutron energy was fixed at 14.7 meV, using a 
horizontally focusing analyzer (HFA). The PG analyzer pieces of HFA are designed to focus 
the final neutron beam, by the pieces forming a spherical surface. Five PG pieces are aligned 
in the horizontal plane and focus the same energy, but different wave vector neutrons. This 
makes the detected intensities almost five times larger than an ordinary flat analyzer. The 
scattering intensities of CEF excitations are almost Q-independent, thus a powder inelastic 
scattering with HFA is an efficient way to obtain sufficient intensities from CEF excitation 
signals. For instance, HFA was used to successfully observe CEF excitations of Ho2Sn207. 4 
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HFA is used not only for CEF excitation experiments, but also used for phonon 5 ' 6 or spin 
excitation measurements under magnetic field 7 studies. 

The conditions of horizontal collimators were 40' - 80' - radial - open. The radial is a radial 
collimator, which divides the neutron beam path by neutron absorbing blades. The sample is 
located on the extension of the blades to optimize the condition of HFA. The Sapphire and 
PG filters were employed to remove higher order contaminations. The energy resolution was 
measured using a vanadium standard sample as AJuo ~ 1.44 meV at the elastic position, 
which is only slightly broader than the calculated value. The 1.1 g pulverized sample was 
sealed in an Al can and mounted in an ILL-type orange cryostat. The constant-Q scans were 
performed at various scattering vectors and temperatures to search the magnetic excitations. 

3. Results 

3.1 Magnetic Structure 

The temperature dependence of the (032) magnetic reflection is shown in Fig. 1. The Neel 
temperature is about 5 K, consistent with the previous reports. 2 ' 3 The magnetic structure, 
the propagation vector and the effective magnetic moments of YbNiSi3 were determined by 
comparing the calculated intensities of several possible models with the measured intensities 
of the magnetic reflections. The magnetic form factor of Yb 3+ was used for these calculations. 8 
The comparison between the calculated and the observed intensities of magnetic reflections 
is shown in Fig. 2. The magnetic reflections were measured at 1.8 K on fifteen reciprocal 
points which are the forbidden positions of the Cmmm space group. The nuclear reflection 
data were used for normalization, where the reflections were measured well above Tr. Figure 3 
represents the revealed magnetic structure of YbNiSi3 . The black arrows denote the directions 
of the magnetic moments. 

The Yb moments are aligned ferromagnetically in the be plane and the ferromagnetic 
layers stacked antiferromagnetically. This magnetic structure is completely different from that 
of the Ce-based compound CeNiGe3, which has the same crystal-structure. 9 The propagation 
vector was determined to be k = (1, 0, 0). The effective magnetic moment of YbNiSi3 was 
determined at 1.1 ± 0.2 ^b- On the other hand, the magnetization measurements at 2 K 
shows that the saturated moment is below 2 u B /f.u. 

Very weak gaussian-shaped reflections were also observed at two forbidden positions above 
Tn- The peak height was about 0.6 cps, while the weakest "not-forbidden" nuclear reflection 
was ~ 1 cps. The reflections at the forbidden positions are assumed to be the results of 
twinning or multiple scattering. If these "forbidden" reflections are caused by twinning, the 
integrated intensities are so weak that the percentage of the twinning should be less than 1%. 
The possibility of multiple scattering has not been confirmed by the experiment. The single 
crystalline YbNiSi3 was so small that it took 15 minutes to measure the (020) reflection, 
a moderately intense fundamental Bragg peak. Although it would be useful to check for 
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the possibility of multiple scattering experimentally, there was not enough time to perform 
additional scans. Those very weak forbidden nuclear reflections were thus neglected when 
resolving the magnetic structure. 




Fig. 1. The temperature dependence of the (032) magnetic reflection. The solid line is a guide to the 
eyes. 




Fig. 2. The calculated (Imag ca i c .) and observed (Imag b s .) intensities of the magnetic reflections. 
The numbers represent the correspond {h k I) reflections. 



3.2 Crystalline Electric Field 

Inelastic neutron scattering measurements were performed at Q = 2.2 and 4.4 A , 
T = 5.5 K (just above Tn) and T = 104 K. Figure 4(a) represents the T-dependence 
for Q = 2.2 A -1 and Fig. 4(b) shows the Q-dependence at T = 5.5 K (just above Tn) 
of the inelastic neutron scattering intensities, which were measured using 5G-PONTA. Three 
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Fig. 3. The magnetic structure of YbNiSi3. The black arrows denote the directions of the magnetic 
moments. The Yb moments align ferromagnetically in the be plane, and the ferromagnetic lay- 
ers stack antiferromagnetically. White, gray and black circles represent Yb, Si and Ni atoms, 
respectively. 

inelastic peaks were observed at AE = 5.5, 12 and 18 meV. The peak heights at 5.5 and 
12 meV are small at high-Q and high-T. Moreover, the Q dependence of the peak heights 
at 5.5 and 12 meV match that of the square of the Yb 3+ form factor reasonably. 8 Thus we 
have confirmed that the peaks observed at 5.5 and 12 meV are indeed magnetic excitations. 
The anomaly at 18 meV is strong at high-Q, while there is no clear temperature dependence. 
Though there might be a contribution of phonon scattering from the Al can 10 or the sample 
itself, more detailed experiments are needed to decide the origin of this anomaly. Quasielastic 
scattering, which is seldom observed in the valence fluctuating systems, were not observed 
nor resolved. At low temperatures, the quasielastic line width has a positive correlation with 
the Kondo temperature. 1 Thus the absence of quasielastic scattering indicates that Kondo 
temperature of this system is considerably small. 

The measured line width of CEF excitation peaks were considerably broader than the 
calculated resolution width. When the peaks at 5.5 and 12 meV were fitted with gaussian 
peaks under the assumption of zero back ground, the line widths (FWHM) were 4.1 ± 0.2 
and 5.5 ± 0.4 meV, respectively. The corresponding calculated resolution widths are 2.2 
and 2.8 meV. Sometimes Ce- and Yb-based compounds show broad CEF excitations. 11-15 
Usually the hybridization of localized 4/ and conduction or valence electrons is employed to 
explain the broadening, while a CEF-phonon coupling was used to explain the broadening for 
YbP04. 15 In addition, for the case of YbP, 16 the broadening was explained by the P-deficiency 
which changes the local symmetry. No indication of CEF-phonon coupling nor distortion of 
Ni or Si were observed in the case of YbNiSi3. Therefore, the hybridization-scenario seems 
most favorable. However, we need more intensive experiment to conclude which scenario can 
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Fig. 4. The results of the inelastic neutron scattering experiments, (a) The temperature dependence. 
The scattering vector was fixed at Q = 2.2 A" 1 , (b) The Q-dependence. The measurements were 
performed at T = 5.5 K (just above Tjj). The solid lines are guides to the eyes. 



4. Discussion 

4-1 Magnetic Structure 

Let us discuss the possibility that the Yb magnetic moments cant toward the a or c 
axis. The related magnetic reflections would not be observed unless the magnetic moments 
cant more than 35° in our experimental condition, mainly due to the smallness of the single 
crystal. However, owing to the representational analysis of the space group Cmmm with 
the propagation vector k = (1, 0, 0), the magnetic structure is determined absolutely in 
the framework of the Landau theory. For a second-order transition, the ordering transition 
can involve only one irreducible representation (IR). The only IR in this case corresponds 
to the model depicted in Fig. 3. In the case of YbNiSi3, the magnetic ordering transition is 
considered to be second order, and no ferromagnetic properties were observed. 2,3 Therefore, 
we have concluded that the magnetic moments do not cant toward any axes from the b axis. 
This representational analysis calculation was performed using the SARA/i-Representational 
Analysis. 17 

4-2 Crystalline Electric Field 

The Stevens' equivalent operator method was used to evaluate the CEF excitation. We 
have assumed a Kramers doublet as the ground state, though it is a subject of a considerable 
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discussion. 2 ' 3 In the orthorhombic symmetry, the CEF hamiltonian is expressed by 6th degree 
polynomial. We used the terms up to 4th degree were used for the analysis of the neutron 
spectrum. That is because it was impossible to reproduce the results of inelastic scattering 
experiments only by 2nd degree terms, while using the terms up to 6th degree made the 
analysis confusing. The CEF hamiltonian which was used to analyze the inelastic scattering 
is 

Hcef = B° 2 Ol + B 2 0\ + B\0\ + B\0\ + BjOt (1) 

where BjS and 0*s are the parameters and the Stevens' operators, respectively. 18 The J = \ 
multiplet of Yb 3+ ion splits into four Kramers doublets under this Hamiltonian. On the 
other hand, we observed only two CEF excitations. Therefore, there might be one more CEF 
excitations that was not observed nor resolved in the present neutron inelastic experiments. 
The derived CEF parameters are listed with the corresponding energy diagram in Fig. 5. 
The sign of Bq was determined negative from the temperature dependence of magnetic sus- 
ceptibility x- 2,19 ' 20 The energy diagram was calculated using the CEF parameters. The black 
arrows denote the excitations that were observed clearly by the inelastic neutron experiments. 
The highest excitation state in Fig. 5 was the result of calculation by the determined CEF 
parameters. 

Though it is hard to determine these parameters absolutely, the physical properties can 
be explained by these determined parameters satisfactorily. The temperature dependence of 
X -1 , which was calculated using those determined CEF parameters in the mean-field approx- 
imation, is presented in Fig. 6 along with the measured value of x~ 1 against log— T? The 
calculated and the measured x~ l against linear— T is plotted in the inset of Fig. 6. In the 
mean- field approximation, is expressed as = \(^EFa — / ^ c " where \ a is a molecular- 
field constant. 21 The value \ ac = \ = — 4.5 mol/emu was used. Although the calculated 
values were slightly different from the observed data in totality, the magnetic susceptibility of 
YbNiSi3 above Tn = 5 K is semi-quantitatively explained by this model. For more complete 
discussion, the 6th degree terms must be included in the CEF hamiltoninan. 

The magnetic contribution of specific heat C mag in an applied magnetic field along the b 
axis was calculated using thus revealed CEF model. The results are shown in Fig. 7. At zero 
magnetic field, a broad peak exists at 20 K, which is considered a Schottky-type contribution. 
This broad peak was not observed in the previous report. 2 ' 3 The reason why the peak has not 
been observed is considered to be the results of the release of entropy by the magnetic order 
transition. Moreover, there is a possibility of the contribution of the lattice specific heat, which 
might make it difficult to resolve the broad peak. When a magnetic field is applied, a sharp 
peak appears at B = 2 T. A possible explanation for this behavior is that the lower lying 
doublet splits by the applied magnetic field. The sharp peak at 2 K becomes broader and moves 
up to higher temperature with increasing the applied field. Above the Neel temperature, these 
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broadenings and movings are qualitatively consistent with the experimental observation. 3 The 
assumption of the Kramers doublet ground state is responsible for these results of calculation. 
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Fig. 5. The level scheme of YbNiSi3. CEF parameters are also presented. The black arrows indicate 
the excitations observed by the neutron inelastic scattering experiments. 
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Fig. 6. The calculated and observed temperature dependence of the magnetic susceptibility against 
log— T. The plots against linear— T is presented in the inset. The solid lines are the calculated 
results. The open circles are the measured values. 2 



Summary 

We studied almost +3 heavy fermion antiferromagnet YbNiSi3 by means of neutron scat- 
tering techniques. We succeeded in determining the magnetic structure and the effective mag- 
netic moment. Moreover, we observed two broad inelastic magnetic peaks, which are attributed 
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Fig. 7. The calculated temperature dependences of the specific heat for various applied magnetic 
fields. Dashed lines are guides to the eyes. 

to CEF excitations. The physical properties of YbNiSi3 were reproduced reasonably well by 
calculation using the CEF level scheme derived from the neutron inelastic scattering. 
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